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f5 ■ Abstract 

O ' We show that strongly monotone systems of ordinary differential equations which have a 

1-g certain translation-invariance property are so that all solutions converge to a unique equilibrium. 
The result may be seen as a dual of a well-known theorem of Mierczynski for systems that satisfy 

G a conservation law. An application to a reaction of interest in biochemistry is provided as an 

>- , illustration. 

^ Keywords: monotone systems, global stability, chemical reaction networks 

_ _ I 

1 Introduction and Motivations 



We recall that a dynamical system is said to be monotone whenever its state space X is endowed 
with a partial order >z and the forward flow preserves the order. In other words, for each ordered 
pair of initial conditions >z ^2, solutions remain ordered: (fti^i) ^ ft{(,2) for all t > 0. See [T3] 
for a discussion and many basic theorems, as well as the recent excellent exposition ^Uj- A 
special and most interesting case is when the partial order is induced by a positivity cone, i.e. a 
closed subset i^T of a Banach space B containing X such that K + K (Z K, K (Z aK for all a > 0, 
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and K U —K = {0}. In this case, one defines a partial order by the rule that ^ ^2 whenever 
^i~^2 £ K. Strict versions of the order are also possible, and particularly useful whenever K has 
non-empty interior: one defines ^1 if ^1 ^ 'C2 and .^1 ^2, and the following even stronger 
notion: ^1 ^ ,^2 if C,i — ^2 ^ int(i^'). A strongly monotone system is one for which the following 
holds: 

6^6 ^ vti^i) » <^t(6) Vt > 0, V6,6 e X. (1) 

A key foundational result is Hirsch's Generic Convergence Theorem ([H |Hl El CHI El), which 
guarantees that, if solutions of such systems are bounded, then, generically, they converge to the 
set of equilibria. Roughly speaking, more complex asymptotic behaviors are possible, but are (if 
they exist at all) confined to a zero-measure set of initial conditions. 

Remarkably, under suitable additional assumptions, generic convergence to equilibria can be 
made global, as is the case if, for instance, the equilibrium is unique [121, sometimes not requiring 
strong monotonicity j^, if the system is cooperative and tridiagonal ^11 or if, there exists 
a positive first-integral for the system, as shown in Mierczinski's paper Our main result 

may be viewed as a dual of the latter result, and applies to strongly monotone systems which 
have the property of translation invariance with respect to a positive vector. Equilibria of such 
systems are never unique. The result is roughly as follows. For systems evolving on Euclidean 
spaces M", we will assume that for some v G int(-ft'), and for all A G R, the following is true: 

ipti^ + Xv) = iptiO + (2) 

for alH G M and all ^ G X. Under strong monotonicity, we show that convergence to equilibria 
is global for a suitable projection of the system. We also show that for competitive systems, i.e. 
systems that are strongly monotone under time-reversal, the same result holds. Statements and 
proofs are in Section |21 

We were originally motivated by proving a global convergence result for certain chemical 
reaction systems which are not necessarily monotone. There has been much interest in recent 
years in estabhshing such global results, see for instance El 11210 CEl III El • In Sectional we 
show how to associate, to any chemical reaction system, a new system of differential equations, 
evolving on a different space (of "reaction coordinates" ) for which our techniques may sometimes 
be applied, and we illustrate with a system of interest in biochemistry. 

In the last section, we make some remarks on extensions and comment on the duality with 
Mierczinski's theorem. 

2 Main Result 

We consider nonlinear dynamical systems of the following form: 

X = fix) (3) 

with states x G X C M", for some closed set X which is the closure of its interior, and some 
locally Lipschitz vector field / : X — > M". For each initial condition G X, we denote by 
^PtiO the corresponding solution, and we assume that V5t(0 is uniquely maximally defined (as an 
element of X) for t G where is an interval in M which contains [0, +00) in its interior. (In 
other words, the system is assumed to be forward -but not necessarily backward- complete.) 



2 



Furthermore, a closed cone K C is given, with non-empty interior, and the corresponding 
non-strict and strict partial orders are considered: y, >~, In particular, we assume that ^ is 
strongly monotone as in and that solutions enjoy the translation invariance property (j2|) for 
some V G int(i^'), which we take, without loss of generality, to have norm one. 

Because of property Q it is natural to assume, and we will do so, that the state space is 
invariant with respect to translation by v, namely: 

xeX^x + XveX vAgM. (4) 

In order to state our main result, we require an additional definition. Given any unit vector 
V, we introduce the linear mapping: 

TT^ : ^ M" : X ^ X - {v'x)v 

(prime indicates transpose), which amounts to subtracting the component along the vector v, that 
is, an orthogonal projection onto f Since {v'x)v = {vv')x, we can also write tt^x = (/ — vv')x. 
Note that ir^v = 0. 

Definition Let ^ G X be given and consider the corresponding solution ft{C,)- We say that <ft{0 
is bounded modulo v if 7T^{ipt{0) is bounded as a function of t, for t >0. 

Notice that we are not asking for precompactness of ^pt{0 (which, in examples, will typically 
fail), but only of its projection. 

Remark Equivalently, the solution ftiO is bounded modulo v if and only if there exists some 
scalar function /3(^, t) : X x [0, +00) M such that ftiO ~ Pi^j ^)'^ is bounded as a function of 
time t. (Recall that X is invariant under translations by v, so this difference is again an element of 
X.) One direction is clear, using /3(^,t) = v'ipt{C). Conversely, suppose that there is any such [3. 
Then: v'^, t)v) = P{^, t)v'v = P{^, t) \v\' = t), so 7r,(/5(e, t)v) = /3(e, t)-{v'm, t)v)v = 0, 
Since X is closed, the assumption is that the closure of {^Pt{0 ~ Pi^^ 1 > 0} is compact. Thus, 
since vr^ is continuous, the same holds for n^^Lpt^^)) = n^^ipt^C,) — /3(^,t)t>). 

We are now ready to state our main result. 

Theorem 1 Consider a forward complete nonlinear system, strongly monotone on X. Let (0) 
enjoy positive translation invariance as in (0) with respect to some vector v G int{K), and so that 
the state space X is closed and invariant with respect to translation by v as in Then, every 
solution which is bounded modulo v is such that TTvi^tiO) converges to an equilibrium. Moreover, 
there is a unique such equilibrium. 

Before addressing the technical steps of the proof, it is useful to provide an infinitesimal 
characterization of translation invariance. This is a routine exercise, but we include a proof for 
ease of reference. 

Lemma 2.1 A system (jS)) enjoys the translation invariance property (j21) with respect to w G M" 
if and only if: 

Xi, X2 G X, Xi — X2 ^ span{f } =^ f{xi) = f{x2) ■ (5) 
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Notice that, for differentiable /, yet another characterization is that v e ker/*(a;) (Jacobian) 
at all states x. 

Proof. If the system is translation-invariant by v, and X2 = Xi + Xv, then (pt{x2) — ^ti^i) — Xv. 
Taking {d/dt)\t=o, we obtain f{xi) — f{x2). We now show the sufficiency of the condition. More 
generally, suppose that L is a linear subspace of R" such that Xi — X2 & L ^ fi^i) — f{^2)', we 
will prove that (ft{x2) — <Pt{xi) is constant if xi — X2 & L. 

We first change coordinates with a linear map T in such a manner that L gets transformed 
into the span of the first i = dimL canonical vectors L = {ci, . . . , e^}. The transformed equations 
are x = f{x), where f{x) = T f{T~^x) and x = Tx. We partition the state as 5; = {y', z')', with 
y of size i, and write the transformed equations in block form: 

y = fi{y,z) 

i = My,z). 

Suppose that two vectors xi and X2 are such that zi — Z2- This means that xi — X2 & L. Then, 
letting Xi := T~^Xi, we have that Xi — X2 G L, and therefore f{xi) = /(.X2) by assumption. 
Thus also f{xi) = Tf{xi) = T f{x2) = /(^2)- In other words, / is independent of y, and the 
transformed equations in block form read: 

y = Mz) 

Z = f2iz). 

Now pick any xi, X2 ^ X such that X1—X2 G L. Then, X1—X2 G L, i.e., Zi = Z2. Let yi{t) and Zi{t) 
denote the components of the solution of the transformed differential equation with respective 
initial conditions Xi, i = 1,2. Then, Zi{t) = Z2(t) for all t > (same initial conditions for the 
second block of variables), which implies that yi{t) = y2{t) for all t. Therefore also Xi(t) = X2{t) 
for all t, and back in the original coordinates we have that {d/ dt){ipt{x2) — <Pt{xi)) = 0, as desired. 

■ 

In order to carry out the proof we first need the following Lemma. 

Lemma 2.2 Let v G int(i^) be given, such that 1^1 = 1. Then, the function: 

V{x) := inf{Q; G R : x ^ av} 
is well defined and Lipschitz for a; G M". 

Proof. We show first that for all x there exists a such that av >z x. We may equivalently check 
that V y x/a for some a 7^ 0. Since x/a ^ as a — > +00, we may conclude that this is the case, 
since, as is well known, {—v,v) := : t; ^ a; ^ —v} is an open neighborhood of the origin, for 
all f ^ (in other words the topology induced by a positivity cone with non-empty interior is 
equivalent to the standard topology in M"). On the other hand, av -< x for all sufficiently small 
a (as q; — > —00, {—x)/{—a) — > 0, so {—x)/{—a) -< v, that is, —x -< —av). Therefore, V{x) 
is well defined. Moreover, since K is closed and the feasible set of a's is bounded from below, 
the infimum is achieved and is actually a minimum, which implies that V{x) is a continuous 
function. We can prove, moreover, that V is Lipschitz, as follows. 
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We first pick an £ > such that ez v for all unit vectors z. (Such an e exists, because 
£2; — i> uniformly on the unit sphere, and (— f , v) is a neighborhood of zero.) Therefore, for each 
two vectors x 7^ applying this observation to z = {x — y), we have that x — y ^ k\x — y\v, 
where k := 1/e, and the same holds if x = y. Now, given any two x, y, we write 

X = x — y + y :< k\x — y\v + V {y)v = {k\x — y\ + V {y))v 

which means that V{x) < k\x — y\ + V{y), and therefore V{x) — V{y) < k\x — y\. Since x and y 
were arbitrary, this proves that V is Lipschitz with constant k. I 

The next Lemma is crucial for proving our main result. 

Lemma 2.3 Let ^1 and ^2 in ^ be arbitrary, and V be defined according to the previous Lemma 
O Then, for alH > it holds that: 

< ^(6-6), (6) 

and the inequality is strict whenever — ^2 ^ span{f }. 

Proof. Let ^1 and ^2 be arbitrary. By definition of V, we have: ^1 ^ ^2+^(^1 ~^2)'W- By translation 
invariance and monotonicity then: (pt{^i) ^ '^t{C,2 + V{C,i,C,2)v) = (pt{^2) + ^(^1,^2)^- It follows 
that V{ipt{^i) — (pt{C,2)) < ^(^1 — ^2), as claimed. In particular, whenever ^1 — .^2 ^ span{t>} 
we have ^1 -< ^2 + V{C,i — C,2)v and therefore, exploiting strong monotonicity: (pt{^i) <^ '^t{C,2 + 
Vi^i - 6)^) = Vti^2) + Vi^i - ^2)v. In particular, then V{ipt{^i) - (^^(6)) < V{^i - 6)- ■ 

Notice that, by the semigroup property for flows. Lemma 12.31 implies that the function t \—>- 
V{Vt{^i) - V^t(6)) is nondecreasing. 

We also prove a result for systems that are strongly monotone in reversed time, meaning 
that for every pair ^1,^2 and every time t < such that (ft{C,i) and <ft{^2) are well-defined the 
following implication holds: 

Corollary 2.4 Let ^1 and ^2 in X be arbitrary, and V be defined according to the previous 
Lemma [2. 2[ Assume that system Q be forward-complete, strongly monotone in reversed time 
over X, and translation invariant with respect to some v G int(-ft'); then, for all t > it holds 
that: 

^(^,(6) - > 1^(6-6), (7) 

and the inequality is strict whenever ^1 — ,^2 ^ span{t>}. □ 

We are now ready to prove the main result. 
Proof of Main Result 

Let ^ G X be such that (pt{0 is bounded modulo v. That is, x{t) := 7iy{ipt{^)) = {I — vv')ipt{C) 
is a bounded function of t. Notice that x satisfies the following differential equation: 

i = (/ - vv')f{^t{0) = (1- vv')f{x) (8) 
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where the last equahty follows by translation invariance. This is a new dynamical system, with 
state space X := {x G f : 3 A G M : x + Af G X}, viz. the projection along w of X onto the 
vector-space v^, and we will denote by ipt the corresponding flow. Notice that vr (we omit the 
subscript v from now on), (pt and ift are related in the following sense: 

Moreover, by translation invariance of X, we have X = X fl f and X = X © span{f }. 

By the above considerations, it makes sense to speak about the tu-limit set uj{x) of solutions 
of ((HI), which by the boundedness assumption, will be a compact, non-empty invariant set. We 
would like to show that uj{x) is a single equilibrium. 

We show uniqueness first. An equilibrium x of (jH)) satisfies that fix) belongs to the span of 
w, let us say fix) = rv. Therefore, the function zit) = x + tf{x) is a solution of the system 
X = f{x), since its derivative satisfies: 

zit) = fix) = fix + irt)v) = fizit)), 

where the second inequality is by (the infinitesimal characterization of) translation invariance. 
Since z{0) = x, we have that (ftix) = x + tf{x) for all t. Assuming that Xi and X2 are two 
distinct equilibria for (jH)), we have that (pt{xi) = Xi + tf{xi) (for i = 1,2). Hence, for all t > 0: 

V{Xi~X2) > V{ipt{Xi) - LftiXi)) 

= Vix,-X2 + [fixi)-fix2)]t) (9) 

By a symmetric argument, however, 

Vix2 - ii) > Vix2 - 5i + [/(52) - fix,)]t) (10) 

which should hold again for all t > 0. It is straightforward, from definition of V{x), that the 
function be increasing with respect to (positive) translations along v. Hence, the inequality in 
(0) implies f{xi) — f{x2) -< 0, while, the second inequality gives /(xi) — /(X2) >- 0. But this is 
clearly a contradiction. 

Let r > be arbitrary; consider the solutions of ^ corresponding to ^ and friO- We claim 
that ftifriO) ~ ftiO is bounded. In fact, denoting by if>t the corresponding projections onto X 
and exploiting Lemma f2.H we obtain: 

^,(^.(0) - vtiO = VrMO) - MO = 1'^^ fivsiO) ds = 1'^^ fiU^m ds (11) 

and so \{pt{(pr{C,)) ~ Vt{0\ — where M is an upper bound on the magnitude of the vector 
field / on a compact set that contains the trajectory Ti^iiptiC))- 

Hence, Viipti^^riyCj) ~ '^tiO) is lower-bounded, and, by virtue of Lemma E31 is decreasing. 
Therefore, it admits a limit V > —00 as t ^ +00. We claim that 

M^riO) - ftiO span{v} . (12) 

Suppose that this claim is false. Then, since, as we just proved, (ptiip^-i^)) — (ftiO is bounded, 
there will be a sequence of times t„ — 00 and an Sq ^ span{T;} such that iptn (v^r(O) ^V'tn (0 ~^ ^o- 
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Moreover, by precompactness of tt o (^i(^), we can pick a subsequence of {tn}, which we denote 
without loss of generahty in the same way, such that vr o ipt^iO ~^ ^o, for some vector xo- 

So the pair [a;o, 5o] belongs to the following set: 

n = {[x,6] : 3tn +00 : iroipt^XO ^ ^ and (fiuX^AO) ' VtniO ^ ^} ■ (13) 
We show next that Q satisfies the following invariance property: 

\/[x,6] en, \ft>0, [^tix),iptix + 5)-iptix)]en. (14) 

Pick any [x, 6] E Q and some sequence {tn} as in the definition of Q, as well as any fixed 
t > 0. From x = lim^^+oo o y^tniO and continuity of the flow, we have: 



^t{x) = lim o = lim n o ip^_^_^^{^) . (15) 

n— >+oo n—>+oo 



Moreover, 



5 = lim (pt„iMO) -Vt„iO 
= lim fr{v>tAO) - ft AO 

n— >+oo 

= hin M^MO) + bVtJO]^) - ^MO) - [v'ftAO]v 

n—f+oo 

where the last equality follows from (^i(7r(^)) = 7i{(pt{0) = VtiO ~ bVilO]"^ applied to t = t„. 
Finally, from the equality (Pr{C + '^'^) = V^r(C) + -^"^ applied to C = <^i„(^(0) and A = v'{pt„{0^ 
this last expression gives that 

6 = lim (/'r(^t„(vr(0)) - <^t„(vr(0) = <^r(5)-x, (16) 

n— >+oo 

that is, a; + 5 = (Pr{x). Therefore: 

<^t{x + 6) - <^t{x) = <^t(<^^(x)) - ^i(x) = lim <^i+r(7ro^i^(^)) -^i(7rov9i„(0). (17) 

71— > + 00 

Now, by translation invariance, we have that: 
and similarly: 
so that: 

so, substituting into (fT7|) . we have: 

ift{x + 6)-(ft{x) = lim (pt+riftAO) - ftiftAO) = lim <ft+tnifriO)-'Pt+t„iO- (18) 

n— >+oo n— >+oo 

Hence, (fT^ follows combining (fT3j) and (fTHj) (using the new sequence {t + 1„}). 
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Recall that V^ipt^ipriO) ~ V^<(0) decreases to its limit V as t ^ oo. On the other hand, for 
any [x,6] G fl, by definition of fl we have that (ft^{ipr{^)) — V^t„(0 S as n oo. Because of 
continuity of V, this implies that V{6) = V. Moreover, by invariance of Q, V {ipt{x + 6) — (pt{x)) = 
V, independently of t. Hence, application of Lemma I2.HI gives 6 G span{f} for any [x,6] G Q. 
This contradicts the assumption that 60 ^ span{f }. Therefore, (I12|l is true. 

Projecting ()12|1 onto the X space shows: 

lim ^tiM-^m ~ 'fMO) = 0. 

t— >+oo 

We next claim that every element of uj{x) is an equilibrium. Indeed, suppose that <^t„(7r(^)) 
p; then, for any r: 

^r{p) = <^r( lim ^t„(7r(0)) = lim (^^ (^t„(7r(0)) = lim ^t„(7r(0) = P- 

Hence, the result follows by uniqueness of the equilibrium for the projected system x = 
il-vv')fix). 

Corollary 2.5 Let a system as in Q be strongly monotone in reverse time and enjoy the 
translation invariance property with respect to some vector v G int(i^'). Then, every solution 
which is bounded modulo v has a projection which converges to an equilibrium. Moreover, there 
is a unique such equilibrium. □ 

Proof. The proof is entirely analogous, once Corollary 12.41 is used in place of Lemma 12.31 I 



3 An Application to Chemical Reactions 

In this section, we show how our result may be applied to conclude global convergence to steady 
states, for certain chemical reactions. A standard form for representing (well-mixed and isother- 
mal) chemical reactions by ordinary differential equations is: 

S = TR{S), (19) 

evolving on the nonnegative orthant ffi>o, where S is an n- vector specifying the concentrations 
of n chemical species, F G M"^™- is the stoichiometry matrix, and R : M>q — M*" is a function 
which provides the vector of reaction rates for any given vector of concentrations. We assume 
that R is locally Lipschitz, so uniqueness of solutions holds, and that the positive orthant M>q 
is invariant, and that it is forward complete: every solution is defined for all t > 0. 

To each system of the form (|TU|) and each fixed vector a G M>o, we associate the following 
system: 

X = f^{x) = R{a + Tx) (20) 

evolving on the state-space 

X, = {x eW^\a + Tx>0} . 
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The ith component Xi of the vector x is sometimes called the "extent" of the ith reaction. We 
will derive conclusions about (fT^ from the study of (pUj) . 

Note that X^- is a closed set which is the closure of its interior (it is, in fact, a polytope), and 
also that X„ is invariant with respect to translation by any v G kerF, because x G means that 
o" + Fx > 0, and therefore also x + Xv E X„ for all A G M, because a + V{x + \v) = a + Fa; > 0. 

As an illustrative example, consider the following set of chemical reactions: 

E+P^C^E+Q 

F + Q^D^F + P, ^ ^ 

which may be thought of as a model of the activation of a protein substrate P by an enzyme E; 
C is an intermediate complex, which dissociates either back into the original components or into 
a product (activated protein) Q and the enzyme. The second reaction transforms Q back into 
P, and is catalyzed by another enzyme (a phosphatase denoted by F). A system of reactions 
of this type is sometimes called a "futile cycle" , and reactions of this type are ubiquitous in cell 
biology. The mass-action kinetics model is obtained as follows. Denoting concentrations with 
the same letters (P, etc) as the species themselves, we introduce the species vector: 

S = {P,Q,E,F,C,Dy 

and these stoichiometry matrix F and vector of reaction rates R{S): 

kiEP - k^iC ' 
k2C 

ksFQ - k^sD ■ 
k^D 

The reaction constants fcj, with i = —1,1,2,3,-3,4, are arbitrary positive real numbers, and 
they quantify the speed of the different reactions. This gives a system (fT^ . Note that, along all 
solutions, one has that 

P{t) + Q{t)+ C{t) + D{t) = constant 

because (1, 1, 0, 0, 1, 1)F = 0. Since the components are nonnegative, this means that, for any 
solution, each of P(t), Q(t), C(t), and D{t) are upper bounded by the constant -P(O) + Q{0) + 
C(0) + D{0). Similarly, we have two more independent conservation laws: 

E{t) + C{t) and F{t) + D{t) 

are also constant along trajectories, so also E and F remain bounded. Therefore, all solutions 
are bounded, and hence, in particular, are defined for all t > 0. The system of equations (jl9|) in 
this example is not monotone, at least with respect to any orthant order. (See [2] for more on 
this example, as well as an alternative way to study it.) We will prove, as a corollary of our main 
theorem, that every solution that starts with E{Q) + C(0) 7^ and F{Q) + -D(O) 7^ converges 
to a steady state, which is unique with respect to the conservation relations. 

Lemma 3.1 The system (j2(J|) is forward complete: every solution is defined for all t > and 
remains in X^j. Furthermore, if it holds that every solution of ()19|) is bounded, then, for every 
solution x{t) of (Uni), Fx()f:) is bounded. 
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Proof. Pick any xq G X^-, and let Sq := a + Txq G M>o- Consider the solution of S{t) of the 
initial value problem S = TR{S), S{0) = So, which is well-defined and satisfies S(t) > for all 
t > 0. Let, for t > 0: 

x{t) := Xo+ [ R{S{T))dT. (22) 
Jo 

Note that x{t) = R{S{t)) for all t. We claim that a; is a solution of a: = fa{x). Since a;(0) = xq 
and X is defined for all t, uniqueness of solutions [fa is locally Lipschitz) will prove the first 
statement of the lemma. To prove the claim, we first introduce the new vector function 

P(t) := a + Tx{t). 

Differentiating with respect to time we obtain that P{t) = Txit) = r{R{S(t))) = S{t) for all 
t > 0. Therefore, P — is constant. Since P(0) = cr + Txq = S{0), it follows that P = S. In 
other words, S satisfies S{t) = cr + Tx(t). Thus, x(t) = R{S(t)) = R{a + Tx(t)) = f^{x(t)), as 
claimed. 

To prove the second statement, we simply remark that, as already proved, for every solution 
X of pO|l . there is a solution S of (|T9|l such that S{t) = a + Tx(t). Therefore, Tx(t) = S(t) — a 
is bounded if Sify is. I 

Note that the futile cycle example discussed earlier satisfies the assumptions of this Lemma. 
We now specialize further, imposing additional conditions also satisfied by the example. 

Lemma 3.2 Suppose that the matrix F has rank exactly — 1, its kernel spanned by some 
positive unit vector v. Let x(t) be a solution of (j2n|l . Then, Vx{f) is bounded if and only if 
7r„x(t) is bounded. 

Proof. Since Ttt^x = r(a; — {v'x)v) = Tx, one implication is clear. Let M be the restriction of F 
to the space f-*- orthogonal to the vector v, i.e. the image of tt^. As Ftt^x = Tx, the images of 
F and M are the same. The map M is one-to-one: suppose that x E is so that if Mx = 0. 
Then, Fx = 0, so x is in the kernel of F, i.e., it is also in the span of v. Thus, x = 0. Let 
be the inverse of M, mapping the image of F into v^. Thus, if a trajectory is such that Fx(t) is 
bounded, then also 

M-^Fx(t) = M-^F7r„x(t) = M-^M'K^x{t) = TT^x{t) 

is bounded. I 

Observe that the spaces X^r are translation invariant with respect to any v as in the statement 
of this Lemma. 

Corollary 3.3 Suppose that: 

1. the matrix F has rank n — 1, with kernel spanned by some positive unit vector; 

2. every solution of (fTUj) is bounded; 

3. (T G M>o is so that the system x = fa{x) is strongly monotone. 
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Then, there is a C = Co- ^ I^>o with the following property: for each p G M>o such that p — cr G 
Image(G), the solution S of (fT^ with S'(O) = p satisfies S{t) — > C as t — > cx). 

Proof. We let the kernel of F be spanned by the positive unit vector v. By Lemmas 13 . II and 13 . 21 
n^xit) is bounded, for every solution of (jSDI)- By Theorem [TJ there is a unique equilibrium ^ of 
the projected system x = [I — vv')f{x) so that every solution x of x = R{a + Fx) is such that 
7Ty{x{t)) C, as t ^ oo. We next show that C = o" + satisfies the requirements. 

Pick p G ]R>o so that p — cr = Fa, a G M™, and let S be the solution of = rR{S) with initial 
condition S{0) = p. Arguing as in the proof of LemmaEHl we have that S{t) = p + Tx{t), where 
X = R{p + Tx), x{0) = 0. 

Introduce the function z{t) = x(t) + a. Then, z = x + = R{p + Fx) = R{a + Tz), with 
z{0) = a. Since a + Fz(0) = cr + Fa = p > 0, it follows that z{0) G X^-, and therefore z{t) is a 
solution of X = R{a + Fx) on X„. Therefore, 7ivz{t) —>■ ^. As x(t) = z{t) — a, this means that 
Hyxit) — > ^ — vr^a. Since for every vector x it holds that Fvr^jX = Fx, applying F to the above 
gives 

Fx(t) = F7r.yx(t) ^ F^ — Fa . 
Therefore, S(t) = p + Tx(t) p + F^ — Fa = cr + F^ = C as t ^ oo. ■ 

In the futile cycle example, we may take f = (1/4, 1/4, 1/4, 1/4)', and consider the following 

set: 

^ = {a = {P,Q,E,F,C,D) eRlol E + C > 0,F + D > 0} . 

The system x = /o-(x) is strongly monotone for cr G S. To see this, we compute the Jacobian of 
i?(cr + Fx(t)) with respect to x: 



(* 


* 





kiE\ 


* 


* 













* 


* 







* 


* / 



where the stars represent strictly positive elements when in the off-diagonals (and strictly negative 
when on the diagonals), and where E, F are the E and F coordinates of cr+Fx, or, more explicitly: 

/* * ki{a3 + (x2 -xi))\ 

* * 

A;3(cr4 + (x4 - X3)) * * 

\0 * * / 

Thus, the system is cooperative (i.e., monotone with respect to the main orthant). It is strongly 
monotone if this matrix is irreducible almost everywhere along trajectories (see e.g. [T^, Sec- 
tion 3.2), which amounts, because fa is a real-analytic function, to asking that cr^ + X2 — xi ^ 
and CT4 + X4 — X3 ^ along any solution. Let us prove now that this is the case, assuming that 
cr G S, that is, that cr^, + ^ and cr4 + erg 7^ 0. Suppose that cr^ + X2 — Xi = 0, so that 
Xi — X2 = and xi — X2 = as. The equations for ^lU^ give: 

xi - X2 = hios + X2- xi)((Ti + X4 - xi) - {k^i + k2){a5Xi - X2) , 

so: 

= -(A;_i + A;2)((T3 + a5) 
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which contradicts as + as 7^ 0. Similarly for (T4 + X4 — X3 = 0. So the system is indeed strongly 
monotone. 

We conclude that every solution of our example with an initial condition in the set S converges 
to an equilibrium. Moreover, there is a unique such equilibrium in each stoichiometry class 
cr + Image(r). 

When initial conditions do not belong to S, one has a standard enzymatic Michaelis-Menten 
type of reaction, and the same conclusion holds. This is very easy to show. (Indeed, take for 
instance the case when E{0) = C(0) = 0. As P = k^D, P{t) is nondecreasing, so (since it is 
upper bounded) we know that P converges. Consider the function y = Q + D. Since P + y is 
constant, y converges, too. Since y has a bounded derivative (it can be expressed in terms of 
bounded variables), and its integral is convergent, it follows ("Barbalat's lemma") that y = —k^D 
converges to zero, so D must converge and therefore, again using that P + Q + D is constant, Q 
converges as well. Finally, since D + F is constant, F converges, too.) 

4 Remarks on Duality and Possible Extensions 

As pointed out in the introduction, our main result stated in Theorem can be seen as a dual 
to Mierczynski's global convergence theorem for strongly cooperative systems with a positive 
first integral, published in [I^j. We discuss this informally in this section. Strictly speaking, 
duality of Q only holds provided that we consider the following special case of Mierczinski's 
Theorem: Consider a system of ordinary differential equations in M", defined by a vector field 
f -.Wl^ M", such that: /(O) = 0, §§- > for all i ^ j, and: 

there exists a vector c G (1^+)" such that d f{x) = for all x G (M+)". (23) 

Then, every solution is bounded and converges to an equilibrium. This is a special case of 
Mierczinski result, which had already appeared in several previous publications, in that linear 
positive first integrals are considered; namely the quantity c'x is preserved along solutions of the 
system. For simplicity, we actually strengthened one of the original assumptions by asking that 
^ > for all i ^ j, rather than a strict monotonicity condition with respect to all off-diagonal 
entries as the Theorem is stated in |13j . 

The duality with Theorem ^ is evident if we express the conditions in terms of the Jacobian 
of the vector field. A linear positive first integral amounts to having a constant left-eigenvector 
relative to the dominant zero eigenvalue for the Jacobian matrix Df{x); in particular, c Df{x) = 
for all X in the state-space. On the other hand, translation invariance by a positive vector v (over 
a given state-space X) can be stated in terms of the Jacobian matrix by asking that Df{x)v = 
for all x G X; i.e., the existence of a constant right-eigenvector relative to the dominant zero 
eigenvalue of the Jacobian matrix Df{x). As a further remark, we note that our main result does 
not need the strict monotonicity condition as stated above, but only asks for strong- monotonicity 
of the resulting flow (this is in fact weaker than assuming strictly positive off-diagonal entries of 
the Jacobian; for instance a much tighter sufficient condition for strong monotonicity of the flow, 
can be formulated by asking that the Jacobian matrix have non-negative off-diagonal entries and 
be irreducible). 
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A more general version of Mierczinski Theorem than stated above does not assume hnearity 
of the first integraL In particular, assumption is replaced by the existence of a function 
H{x) -.Wl^R, such that DH{x)-f{x) = and DH{x) e M+ for all x G M+. This condition does 
not allow an elegant interpretation in terms of Jacobians of Df{x), but nevertheless, one may 
state a nonlinear dual of the Theorem would provided that we understand translation invariance 
in the following more general sense. Let us say that a flow is invariant with respect to translation 
by a strictly increasing flow (p if for all tx.t^ in R the following holds: 

<^ti(^t2(a;o)) = <^t2(<^ti(a;o)), 

and moreover, for each Xi,X2 € X there exists t G M so that X2 ^ ipt{xi). This property 
generalizes our previous concept: translation invariance with respect to a constant vector v is 
exactly the property of invariance with respect to translation by the increasing flow if induced 
by the system of differential equations x = v. Invariance with respect to non-trivial general flows 
as in this definition is not easy to check in concrete examples, however, at least in principle, an 
infinitesimal characterization of the property is as follows. Let f{x) : X — and v{x) : X 
be vector- fields. The flow induced by the system x = f{x) commutes with respect to the 
strictly increasing flow induced by x = v{x) if and only if: 

Df\x)v{x) = Dv{x)f{x) . 

Moreover, if there exists a compact set P C mt{K) so that v{x) G P for all x G X, then the flow 
induced by v{x) is strictly increasing (meaning that its solutions are such with respect to t) and 
for any Xi and X2 in X there exists t G M so that X2 ^ ipt{xi). 

Accordingly we have to redefine the notion of boundedness modulo translation by (f by asking 
that solutions are bounded if there exists M > such that for all xq G X and all t G M, there 
exists r with the property that \(fr{ft{xo))\ < M. While this definition is rather natural, there 
is not, however, a natural counter-part to the space X = X fl v^. Hence, we may as well let 
X be defined as a quotient space of X/ ~ under the equivalence relation xi ~ X2 if and only 
if <ft{xi) = X2 for some t G M. This definition of X and the commutativity of ip and allow 
us to define a flow on equivalence classes of [x] of X in the natural way: </3<([x]) := [y9t(x)]. 
Boundedness of a solution in the space X is equivalent to boundedness modulo translation given 
above. Our main result would then be translated into the following statement in the current 
set-up: Consider a forward complete, strongly monotone nonlinear system (0) with translation 
invariance with respect to a strictly increasing flow. Then, every solution which is hounded is 
such that ^t{\x\) admits a limit as t ^ +00. 
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